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1. The model 

T— I ■ 

I We present a detailed, physically grounded, model for the early co- 

■ evolution of spheroidal galaxies and of active nuclei at their centers. 
T^lj- , The model is based on very simple recipes, that can be easily imple- 
0^ ' mented. The main components and the transfer processes accounted for 

■ are depicted in Fig. 1, while we defer to Granato et al. (2004) for a full 
I description. In summary, we start from the diffuse gas within the dark 

matter halo falling down into the star forming regions at a rate ruled by 
Q^i the dynamic and the cooling times. Part of this gas condenses into stars, 

'I at a rate again controlled by the local dynamic and cooling times. But 

^ • the gas also feels the feedback from supernovae and from active nuclei, 

, heating it and possibly expelling it from the potential well. Also the ra- 

' diation drag on the cold gas decreases its angular momentum, causing an 

J> . inflow into a reservoir around the central black hole. Viscous drag then 

^ I causes the gas to flow from the reservoir into the black hole, increasing 

^ • its mass and powering the nuclear activity. Among the most novel point 

of our work with respect to other semi-analytic approaches, we point out 
the treatment of the evolution of QSO activity in the forming spheroids, 
and of the resulting effects on galaxy evolution. 

2. Results 

In the shallower potential wells (corresponding to lower halo masses 
and, for given mass, to lower virialization redshifts), the supernova heat- 
ing is increasingly effective in slowing down the star formation and in 
driving gas outflows, resulting in an increase of star/dark- matter ratio 
with increasing halo mass. As a consequence, the star formation is faster 
within the most massive halos, and the more so if they virialize at sub- 
stantial redshifts. Thus, in keeping with the proposition by Granato et 
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Figure 1. Scheme of the baryonic components included in the model (boxes), and 
of the corresponding mass transfer processes (arrows). 



al. (2001), physical processes acting on baryons effectively reverse the 
order of formation of galaxies compared to that of dark-matter halos. 

A higher star-formation rate also implies a higher radiation drag, 
resulting in a faster loss of angular momentum of the gas (Umemura 
2001) and, consequently, in a faster inflow towards the central black- 
hole. In turn, the kinetic energy carried by outflows driven by active 
nuclei through line acceleration is proportional to -^^Qgo ^'^^ this mech- 
anism can inject in the interstellar medium a sufficient amount of energy 
to unbind it. The time required to sweep out the interstellar medium, 
thus halting both the star formation and the black-hole growth, is again 
shorter for larger halos. For the most massive galaxies (Myir ^ 10^^ Mq) 
virializing at 3 ;$ z^ir ;$ 6, this time is < 1 Gyr, so that the bulk of the 
star-formation may be completed before type la supernovae can signifi- 
cantly increase the Fe abundance of the interstellar medium; this process 
can then account for the a-enhancement seen in the largest galaxies. 

The interplay between star formation and nuclear activity determines 
the relationship between the black-hole mass and the mass, or velocity 
dispersion, of the host galaxy, as well as the black-hole mass function. As 
illustrated by Figs. 2 and 3, the model predictions are in excellent agree- 
ment with the observational data. A specific prediction of the model is 
a substantial steepening of the M^y{-(t relation for a ^ 150 km s^^: the 
mass of the BH associated to less massive halos is lower than expected 
from an extrapolation from higher masses, because of the combined ef- 
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Figure 2. Predicted relationship between black-hole mass and line-of-sight velocity 
dispersion of the host galaxy for different virialization redshifts. 
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Figure 3. Predicted local black-hole mass function compared with the recent esti- 
mate by Shankar at al. (in preparation). The total mass density in BHs they derive 
is pBH ~ 4 X IO^MqMpc"^, 25% less than our model. The decline of the model at 
low Mbh is due to having considered only objects with Mvir > 2.5 x lO^^M©. 

feet of SN heating, which is increasingly effective with decreasing galaxy 
mass in hindering the gas inflow towards the central BH, and of the 
decreased radiation drag. 

The ratio between QSO and SNae feedback is an increasing function of 
the mass. At low mass the cumulative effect of QSO is almost negligible 
with respect to that of SNae, but it becomes dominant at intermedi- 
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ate mass (typically by a factor of a few) and high mass (by a factor 
^ 10). Note that the QSO effect usually increases exponentially with 
time, while that of SNae increases more slowly. Thus the instantaneous 
QSO effect becomes dominant, if ever, only a few e-folding times before 
the maximum of QSO activity. 

Coupling the model with GRASIL (Silva et al. 1998), the code com- 
puting in a self-consistent way the chemical and spcctrophotometric evo- 
lution of galaxies over a very wide wavelength interval, we have obtained 
predictions for the sub-mm counts and the corresponding redshift dis- 
tributions as well as for the redshift distributions of sources detected 
by deep K-band surveys, which proved to be extremely challenging for 
all the current semi-analytic models. The results, shown by Fig. 4, are 
again very encouraging. 

The cosmic history of BH accretion, is in keeping with the common no- 
tion that quasar activity peaks at redshift 2-3. Also, the formation rate 
of spheroids derived in this paper is consistent with that we computed 
in Granato et al. 2001 from a deconvolution of the high-z luminosity 
function of QSO. 




Figure 4- Left: predicted 850 /xm cxtragalactic counts compared with observed 
SCUBA counts. Right: predicted redshift distribution of galaxies brighter than K = 
20 compared with the results of the K20 survey 
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